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ABSTRACT 

We have discovered an ultrasoft X-ray transient source, 2XMMi J184725. 1-631724, which was de- 
tected serendipitously in two XMM-Newton observations in the direction of the center of the galaxy 
IC 4765-f01-1504 at a redshift of 0.0353. These two observations were separated by 211 days, with 
the 0.2-10 keV absorbed flux increasing by a factor of about 9. Their spectra are best described 
by a model dominated by a thermal disk or a single-temperature blackbody component (contributing 
>80% of the flux) plus a weak power-law component. The thermal emission has a temperature of a few 
tens of eV, and the weak power-law component has a photon index of ^3.5. Similar to the black hole 
X-ray binaries in the thermal state, our source exhibits an accretion disk whose luminosity appears 
to follow the L oc T 4 relation. This would indicate that the black hole mass is about 10 5 -10 6 M Q 
using the best-fitting inner disk radius. Both XMM-Newton observations show variability of about 
21% on timescales of hours, which can be explained as due to fast variations in the mass accretion 
rate. The source was not detected by ROSAT in an observation in 1992, indicating a variability factor 
of >64 over longer timescales. The source was not detected again in X-rays in a Swift observation in 
2011 February, implying a flux decrease by a factor of >12 since the last XMM-Newton observation. 
The transient nature, in addition to the extreme softness of the X-ray spectra and the inactivity of 
the galaxy implied by the lack of strong optical emission lines, makes it a candidate tidal disruption 
event. If this is the case, the first XMM-Newton observation would have been in the rising phase, and 
the second one in the decay phase. 

Subject headings: accretion, accretion disks — galaxies: individual:2XMMi J184725. 1-631724 — galax- 
ies:nuclei — X-rays: galaxies. 



1. INTRODUCTION 

It is generally believed that supermassive black holes 
(SMBHs) reside in many galaxies, but only a fraction 
of them exhibit a ctive galactic nuclei (AGN), w hile oth- 
ers are dormant (|Kormendv fc Richstond Il995| ). Tidal 
disruption events can provid e a unique way to find and 
study the dormant SMBHs (|Reeslll988ft . Such an event 
occurs when a star approaches a SM BH and is tidally 
disrupted and s ubsequently accreted ()Lidskii fc OzernoH 
[I^79l:fReei[T988h . The mass of the SMBH should be <10 8 
Mq for such events to occur outside the event horizon for 
solar-type stars. Tidal disruption events are expected to 
be transient and rare, wi th the avera ge occurrence rate of 
~ 10~ 4 yr _1 per galaxy (|Reeslll99"ol ). They are predicted 
to have a fast rise, with a timescale of half a year, and 
the decay can last on the order of months to years, with 
the luminosity decayin g as L oc i~ 5 / 3 ()Lidskii fc Ozernol 
[1971 |R^[1^M[1990] ). The peak of the flare is expected 
to reach the Eddington luminosity and be dominated by 
thermal UV or X-ray emission. 

A few tidal disruption event candidates were 
found from the ROSAT All-Sky Survey, such as 
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RX J 1624.9+7554 RX J1242.6-1119 and NGC 
5905 dGrupe et all [l999t E omossa fc Greinerl |l999; 
iKomossa &: Badelll999D . Their host galaxies were con- 
firmed to be inactive or only weakly active (NGC 5905) 
using Hubble Space Telescope spectroscopy ()Gezari et al.l 
2003;. They had peak soft X-ray luminosities up to 
~10 44 erg s _1 and showed ultrasoft X-ray spe ctra with 
blackbody temperatures of ~0. 04-0.1 keV (|Komossal 
l2lMl2TjQ8l ). NGC 5905 is the best observed candidate 
and is the first one found to follow approximately the 
L oc i~ 5 / 3 evolution. Some candidates wer e detected re- 
cently from the XMM-Newton Slew Survey (|Esquei et al.l 
I2007ll200l) in the X -rays and from the Gal axy Evolution 
Explorer D e ep Im aging Survey in the UV ()Gezari et al.l 
120061 120081 I2009D . Very recently, the transient source 
Swift J164449. 3+5734 51 was suspected to be due to a 
tidal disruption event (jBloom et al.ll201lHBurrows etafl 
2011). In contrast with the candidates above and the 
theoretical prediction, this sourc e is hard in X-rays , 
with a photon index around 1.8 (|Burrows eTalll2l)ll . 
Its tidal disruption event explanation still needs to be 
confirmed by future long-term monitoring. 

We are carrying out a project of classifying a sam- 
ple of sources in the Second XMM-Newton Ser endipi- 
tous Source (2XMM) Catalog (jWatson et al.ll2009[ ). Here 
we report on the discovery of an ultrasoft X-ray tran- 
sient source, 2XMMi J184725. 1-631724, whose position 
is RA=18:47:25.16, Dec=-63:17:24.96 (J2000) from the 
2XMM catalog, with a 1-ct error of 0726. It is in the 
direction of the cent er of the galaxy IC 4765-f01-1504 
(jCarrasco et al.|[2006f ). This source has negligible emis- 
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sion above 2 keV. We describe the multi-wavelength ob- 
servations of the source and the data reduction in Sec- 
tion [2] In Section |3l we first give the multi- wavelength 
detections of the source, followed by presentations of its 
detailed X-ray spectral and timing properties. We dis- 
cuss its possible nature in Section |4] and draw our con- 
clusions in Section [SJ 

2. DATA ANALYSIS 

2.1. XMM-Newton Observations 

2XMMi J184725. 1-631724 was observed twice by 
XMM-Newton (Tabled)), on 2006 September 7 and 2007 
April 16. These two observations of this source will 
be referred to hereafter as XMM1 and XMM2, respec- 
tively. The source was detected in all the three Euro- 
pean Photon Imaging Cameras in the imaging mode, i.e., 
pn, M OS1, and MOS2 (IJansen et a l.ll2001b IStriider et al.l 
I200H iTurner et all l2001lh in both observations. The 
source was also d etected by the Optical Monitor (OM; 
iMason et aTll200lD in XMM1, but it was not in the FOV 
of the OM in XMM2. In XMM1, the two UV filters 
UVW1 and UVM2 were used, and we obtained the source 
detection information directly from the pipeline prod- 
ucts. 

We used SAS 10.0.0 and the calibration files of 2010 
November for reprocessing the X-ray event files and 
follow-up analysis. The data in strong background flare 
intervals, mostly at the end of the XMM2 observation in 
the pn camera, are excluded following the SAS thread for 
the filtering against high backgrounds. The final expo- 
sures used are given in Table [TJ We extracted the source 
spectra of the pn, MOS1, and MOS2 cameras from a cir- 
cular region centered on the source using 15" and 35" 
radii for XMM1 and XMM2, respectively. A smaller ra- 
dius was used for XMM1 because the source was fainter 
and near the CCD gap. The background spectrum was 
extracted from a large circular region with a radius of 
100" near the source in each camera. The event selec- 
tion criteria fol lowed the default va lues in the pipeline 
(see Table 5 in iWatson et al.1 ()2009f )). We rebinned the 
spectra to have at least 20 counts in each bin so as to 
adopt the \ 2 statistic for the spectral fits. 

We also extracted light curves from the pn camera, 
which has a larger effective area and a higher timing res- 
olution than the MOS cameras, using the same aper- 
tures as those for spectral extraction. We first extracted 
background-subtracted light curves with a bin size of 
250 s, using the SAS task epiclccorr to apply relative 
corrections. To create the power density spectra (PDS), 
we also extracted light curves from the source region us- 
ing the frame time as the bin size, which is 199.1 ms for 
XMM1 (using the extended-full-frame mode) and 73.4 
ms for XMM2 (using the full- frame mode). Considering 
that the source is very soft and the background domi- 
nates above 2 keV, all light curves were extracted in the 
energy range 0.2-2.0 keV. We calculated the PD S us- 
ing a similar procedure as, e.g.. iGoad et a l. (2006). The 
XMM1 199.1 ms and XMM2 73.4 ms pn light curves were 
split into segments each with 32768 and 65536 data bins, 
respectively, resulting in four segments for XMM1 and 
five for XMM2. The PDS was calculated for each seg- 
ment, and all PDS for each light curve were merged and 
averaged by binning in frequency using a logarithmic fac- 



tor of 1.1, under the condition that each bin contains at 
least 20 individual PDS measurements. The errors were 
calculated from the sample standard deviation of PDS 
measurements in each bin. 

2.2. ROSAT and Swift Observations 

Our source was not detected in the ROSAT All-Sky 
Survey in 1990, which had a d etection limit of .1-2.4 
keVflux 5xl0" 13 erg s" 1 cm" 2 (jVoges et al.l ll999h Our 
source was in the FOV of one ROSAT PSPC pointed ob- 
servation (the sequence number 800256, 1992 October, 
~11 ks), at an off-axis angle of ~2.6'. It was not detected 
either and was (thus) not listed in the WGA catalog of 
the R OS A T point sources (jWhite et al.lll994f ). We calcu- 
lated the confidence interval of the s ource detection us - 
ing Bayesian statistics as described in Kraft e t al.1 (|1991| ). 
Circular source and background regions with radii of 40" 
and 2' respectively were used. The corresponding (an- 
cillary plus photon redistribution) response matrix was 
generated and used to convert the count rates to the 
fluxes. 

At our requ est, the Swift Gam ma Ray Burst Ex- 
plorer mission (IGehrels et al.|[2004T) observed the field of 
2XMMi J184725.1-631724 on 2011 February 23 for a to- 
tal of 5 ks (obs ervation ID 00031930 001). The X-ray 
telescope (XRT; iBurrows et al.l I2005D wa s operated in 
Photon Counting mode ()Hill et al.l I2004D . X-ray data 
were reduced with the task xrtpipeline version 0.12.1. 
We found an enhanced count rate at the position of our 
source, but it is very weak. We also calculated the con- 
fidence interval of the detection. Radii of 23"5 and 235" 
were used for the circular source and background regions, 
respectively. The corresponding response matrix was 
generated using the calibration fil es of 2011 February . 
The UV-Optical Telescope (UVOT; IRoming et al|[2005l) 
was operated using the UVW1 filter for 5 ks. The mag- 
nitude and flux were measured with the task uvotsource 
version 3 based on the mos t recen t UV OT calibration 
as des cribed in iPoole et"aLl (|2008l ) and Brccvcl d et al.l 
(2010). Circular source and background regions with 
radii of 5" and 20", respectively, were used. 

2.3. Optical Observations 

Our source is in the direction of the cen ter of the galaxy 
IC 4765-f01-1504 (jCarrasco et al.|[2006h . This galaxy is 
located in the background of the rich group of galaxies 
IC 4765 (also known as Abell S0805, z=0.01497). It was 
imaged with the 1.3 m Warsaw telescope at Las Cam- 
panas Observatory in Chile through the standard John- 
son V and Cousins I filters in 1998. We us ed the V- and 
I-filter images from iCarrasco et al.l (|2006l ) to derive the 
main photometric p arameters of the g alaxy with a Sersic 
model in GALFIT (|Peng et al.l l2010t The images have 
a FWHM of the P SF of about 1"2. 

ICarrasco et al.l (|2006| ) also obtained an optical spec- 
trum of the galaxy on 1999 June 19 with the Wide Field 
CCD camera mounted on the 2.5 m Du Pont Telescope 
at the Las Campanas Observatory in Chile, but it has 
poor quality. We obtained a new longslit spectrum of 
this gala xy with the Gemin i Multi-Object Spectrograph 
(CMOS, IHook et al.1 [200l at the Gemini South Tele- 
scope in the queue mode. The observation was made on 
the night of 2011 March 19 (UT) during bright time (il- 
lumination fraction 0.99), under photometric conditions 
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TABLE 1 
XMM-Newton Observation Log 



Observation ID Date off-axis angles (arcmin) Duration Exposurc(ks) Filter 

pn/MOSl/MOS2 ' (ks) pn/MOSl/MOS2 



0405550401(XMM1) 2006-09-06.98 3.9/3.0/3.7 28.0 19.5/27.6/27.6 medium 

0405380501(XMM2) 2007-04-16.31 9.0/8.5/9.4 34.7 20.5/32.2/31.9 thinl 



and rsj\" seeing. The 400 lines/mm ruling density grating 
(R400) centered at 5500A was used, to minimize the ef- 
fect of moon illumination. The slit width was set to 1". 
A total of four exposures of 900 s each were obtained. 
Small offsets in the spectral direction (50A) towards the 
blue and the red were applied between exposures to al- 
low for the gaps between CCDs and to avoid any loss 
of important lines present in the spectra. Spectroscopic 
flats and comparison lamp (CuAr) spectra were taken 
after each science exposure. In addition, the spectropho- 
tometric standard star LTT 7379 was observed at the 
end of the night to flux calibrate the science spectrum. 

The observations were processed with the Gemini 
IRAF package version 1.9 in IRAF. All science expo- 
sures, comparison lamps and spectroscopic flats were bias 
subtracted and trimmed. Spectroscopic flats were pro- 
cessed by removing the calibration unit plus CMOS spec- 
tral response and the calibration unit uneven illumina- 
tion, normalizing and leaving only the pixel-to-pixel vari- 
ations and the fringing. The resulting two-dimensional 
spectra were then wavelength calibrated, corrected by 
S-shape distortions, sky-subtracted, extracted to a one- 
dimensional format using a fixed aperture of 7'.'8 in di- 
ameter, and then average combined. The final spectrum 
has a resolution of ^8.8 A (FWHM) and a dispersion 
of ^1.36 A pixel -1 , covering a wavelength interval of 
^4000- 7600 A. The signal-to-noise ratio is about 40 at 
5500 A. 

We measured the redshift of the galaxy with two meth- 
ods. In the first method, we cross-correlated the spec- 
trum with a high signal-to-noise template using the fx- 
cor routine in the IRAF RV p ackage. The error wa s 
estimated using the R statistic of iTonrv k, Davisl (| 19791 ): 
a v = (3/8)(w/(l + i?)), where w is the FWHM of the cor- 
relation peak and R is the ratio of the correlation peak 
height to the amplitude of the antisymmetric noise. In 
the second method, we identified the most prominent ab- 
sorption lines (as no clear emission lines were detected) 
in the spectrum and derived the redshift by employing a 
line-by-line Gaussian fit using the rvidline routine in the 
IRAF RV package. 

3. RESULTS 

3.1. The Source and the Multi-wavelength Observations 

2XMMi J184725. 1-631724 was detected in X-rays in 
both XMM1 and XMM2. We see no clear X-ray emis- 
sion of our source from the R OSAT observation in 1992 
October and the Swift observation in 2011 February, in- 
dicating variability factors of >64 and >12, respectively, 
compared with XMM2 (the 0.2-10 keV absorbed flux; 
see Section l3~4l) . Here we concentrate on observations in 
other wavelengths. 

The source position from the 2XMM catalog has been 
astrometrically corr ected by matching with the optical 
catalog USNO B1.0 (|Monet et al.ll2003l) . For XMM1 and 
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Fig. 1. — The optical image of the galaxy in the V band. The 
pixel size is fy/414. The green plus marks the central position of the 
galaxy obtained from the fits to its V- and I-band profiles using 
a Sersic model, and it is at RA=18:47:25.14, Dec=-63:17:24.77 
(J2000). The red circle is centered at the X-ray position, with the 
radius corresponding to the 3-cr error. 

XMM2, a very small correction was invoked (a fraction of 
an arcsec) . We compared the corrected positions of ten of 
the brightest X-ray sources in each observation with the 
positions of optical counterparts from the USNO B1.0 
catalog and found that most of the offsets between the 
matches are less than 0.5", indicating successful astro- 
metric corrections. 

Table [5] gives the detection of a UV source in UVW1 
and UVW2 from the XMM1 OM at a position within 
the 2-a positional error from 2XMMi J184725. 1-631724 
and is probably its UV counterpart. There is also a UV 
source detected near our source in UVW1 from the Swift 
UVOT. After applying the astrometric correction using 
the USNO-B1.0 catalog with the XMM-Newton SAS task 
eposcorr, we obtain its position of RA=18:47:25.14 and 
Dec=-63:17:25.04 (J2000), within the 1-a error from the 
X-ray position. Its magnitude and flux are 18.67±0.04 
and (1.34±0.05)xl0~ 16 erg s -1 cm -2 A -1 , respectively. 
We note that the UV filter set in the UVOT is different 
from that of the XMM-Newton OM, and the above values 
should not be directly c ompared with the OM measure- 
ments in XMM1 above. IGrupe et al.l (|2008aD measured 
an offset between the magnitudes from the two instru- 
ments of W^Iom — W^Iuvot = 0.78 by comparing several 
field stars in the images of the AGN Mkn 335. With this 
taken into account, there seems to be little variability in 
the UV between the two epochs. 
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TABLE 2 

The counterpart candidates in UV, optical, and IR 



UV (XMM1 OM) 


Optical (USNO B1.0) 


IR (2MASS PSC) 


r UVW1 UVM2 
(arcscc) (AB mag/flux(lCr 16 erg s _1 cm~ 2 A -1 )) 


r B2 R2 I 

(arcsec) (mag) 


r J H K 

(arcsec) (mag) 


0.52 19.65±0.10/1.78±0.17 19.83±0.19/2.38±0.41 


0.16 16.1 15.33 15.40 


0.31 15.29±0.07 14.59±0.09 14.30±0.08 



Note. — The r columns are the offsets of the counterparts from 2XMMi J1S4725. 1-631724. The magnitudes/fluxes are not 
corrected for the Galactic reddening. 
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Fig. 2. — The smoothed spectrum of the galaxy IC 4765-f01-1504 from the Gemini South Telescope, with the most important absorption 
lines identified. The drop at 6870 A is due to the atmosphere OH absorption. 



The optical and IR counterpart candidates of the 
source from the USNO B1.0 and 2MASS Point Source 
Catalogs are given in Table [2] The optical counter- 
part, the galaxy IC 4765- f01-1504, is shown in Fig- 
ure Q] (jCarrasco et alJl2006f ). Our fits of the V- and I- 
band images using a Sersic model give integrated mag- 
nitudes of 16.99±0.01 and 15.62±0.02, effective radii of 
2'/52±0'.'04 and 2'.'94±0"09, Sersic indices of 3.54±0.06 
and 4.14±0.13, and apparent axis ratios of 0.238±0.003 
and 0.271±0.003, respectively. This galaxy may be an 
elliptical galaxy, which typically has a Sersic index of 4. 
The axis ratios above would imply a high inclination of 
this galaxy if its intrinsic ellipticity is low. 

Figure [5] shows the spectrum of IC 4765-f01-1504 from 
the Gemini South Telescope. No clear emission lines 
were detected, supporting the identification as an ellipti- 
cal galaxy. A redshift of z=0.0353±0.0001 was obtained 
from both the cross-correlation method, with R — 12.80, 
and the absorption line fit method, indicating a perfect 
agreement between them. This redshift disagrees wit h 
the value of 0.0869 obtained by iCarrasco et all (|2006l) . 
which used the cross-correlation method (as there were 
no significant emission lines detected). Considering that 
our new spectrum has much better quality, we adopt this 
new redshift. We measured a 3-cr upper limit of [OIII] 
5007 Aof0.9xl0- 15 erg s 1 cm 2 . Assuming a flat uni- 
verse with the Hubble constant H =73 km s _1 Mpc -1 
and the matter density f2M=0.27, this redshift corre- 



sponds to a comoving radial distance of 143.9 Mpc and 
a luminosity distance of 149.0 Mpc, which will be used 
in this paper. The absolute V and K magnitudes of this 
galaxy are -19.2 and -21. 6, respectively, after the Galac- 
tic extinction correction (jSchlegel et al.lll998l ). Based on 
the BH mass vs. bulge lum inosi ty relations fromlGrahaml 
(|2l^ . lLauer et all pOOl and lMarconi fc Hu nt ( 2003), 
the above magnitudes imply the mass of the SMBH in 
IC 4765-f01-1504 to be about 10 7 and 10 6 M if the 
bulge/total luminosity ratio is 1 or 0.1, respectively. Be- 
cause the sample of galaxies in the above studies were 
generally brighter than IC 4765-f01-1504 and these rela- 
tions have large intrinsic scattering, these mass estimates 
might have an uncertainty as large as one order of mag- 
nitude. 

3.2. X-ray Spectral Modeling 

We fitted the spectra of 2XMMi J184725. 1-631724 
from both XMM1 and XMM2 using various spectral 
models. We jointly fitted the spectra from all three cam- 
eras, i.e., pn, MOS1, and MOS2, and their relative nor- 
malizations were left free. We only report the normaliza- 
tion results corresponding to the pn camera. MOS1 and 
MOS2 differ by about 10% (the largest one -20%), less 
than the error bars. We fitted the spectra in the 0.2-10 
keV energy band. 

We were unsure the X-ray emission mechanism of 
our source. Thus we first tested the common single- 
component models to see whether any of them can de- 
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TABLE 3 
Spectral modeling results 



Model 


Obs 


JVh 

(10 20 cm" 2 ) 


^MCD/BB 

(eV) 


^MCD/BB 

(10 4 ) 


TpL/SIMPL 


Npl//sc 
(10" 5 /%) 




/mcd/bb 

(%) 


(10~ 12 org 


■Funabs 
s- 1 cm" 2 ) 


L 

(10 43 org s" 1 ) 


MCD+PL 


XMMl 
XMM2 


8.6±0.5 


65.8±5.0 
93.1±2.2 


1 5fi+ 121 
1 45+ 034 


3.72±0.62 
3.27±0.70 


1.79±0.32 
3.75±1.39 


0.92(115) 
1.20(284) 


96.4+i;" 


0.22±0.01 
1.93±0.03 


1.76±0.19 

io.i8±;-»? 


1.70±0.44 
6.38±0.66 


SIMPL(MCD) 


XMMl 
XMM2 


8.5±0.5 


65.1±5.5 
92.9±2.4 


1 7Q-t-l't>0 
1 - ' J -0.73 

1.49±0.32 


3.71±0.63 
3.33±0.78 


o got'* 


0.92(115) 
1.20(284) 


90.0 + 7 I 
97.7l°- 9 2 


0.22±0.01 
1.92±0.03 


1. < O_ 24 

10.21±0.89 


1.81±0.57 
6.40±0.72 


BB+PL 


XMMl 
XMM2 


7.4±0.7 


57.8±3.9 
78.2±1.6 


1 63 +1 ' la 
2.37±0.41 


3.71±0.59 
3.89±0.61 


1.76±0.28 0.92(115) 
5.54±1.46 1.18(284) 


78.2+- 5 2 
86.0+I 4 6 4 


0.22±0.01 
1.92±0.03 


1.23±0.11 

7 - 59 io.51 


0.59±0.14 
2.79±0.39 


SIMPL(BB) 


XMMl 
XMM2 


7.2±0.5 


57.7±4.1 
77.5±2.2 


l Ki+ la 

1 - 81 -0.65 

9 ^S+ 0/rl 

2 - t>8 -0.45 


3.71±0.59 
4.25±0.83 


5.48i 4 .- 

q oc + 3.56 
°' 00 — 1.63 


0.92(115) 
1.18(284) 


94.2t 2 ; 4 


0.22±0.01 
1.92±0.03 


1.18±0.18 
7.2±0.68 


0.59±0.15 
2.70±0.31 



Note. — The column /mcd/bb refers to the unabsorbed flux fraction of the MCD/BB component in the 0.2—10 keV energy band (for the 
models SIMPL(MCD) and SIMPL(BB), it refers to the unscattcrcd part). -F a bs and .Funabs are the total absorbed and unabsorbed fluxes in the 
0.2—10 keV energy band, respectively. The luminosity L was calculated using the unabsorbed bolomctric flux of each spectral component (the PL 
component was integrated down to 0.2 keV). All errors are at a 90%-confidcnce level. 
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Fig. 3. — The unfolded spectra and the fit residuals using the model MCD+PL. For clarity, only the pn spectra are shown for the unfolded 
spectra. The dotted, dashed, and solid lines are for the MCD and PL components and the total model, respectively. The residuals are 
shown for all three cameras (black/red/green for pn/MOSl/MOS2, respectively). 



scribe our X-ray spectra well: a single temperature black- 
body (BB), a multi-color disk (MCD), a PL, a broken 
PL, a cut-off PL, an APEC thermal plasma model, and 
a thermal bremsstrahlung spectrum. They are models 
bbodyrad, diskbb, powerlaw, bknpower, cutoffpl, APEC, 
and bremss in XSPEC, respectively. All models include 
the absorption described by the WABS model in XSPEC; 
our results change little with alternative absorption mod- 
els such as PHABS or TBABS in XSPEC. All these sim- 
ple models fail to describe one or both of the XMM- 
Newton spectra, with residuals above 1 keV typically 
seen. For indication, we report the PL index Tpl of the 
fits using the PL model. We obtained F PL = 5.86±0.27 
for XMMl (x^(V)=1.13(117)) and 6.88±0.09 for XMM2 
((X^0)=2.25(286)). The lower X 2 value for XMMl to 
some degree is due to poorer data. 

We next attempted to fit the spectra with the double- 
component models MCD+PL and BB+PL, finding that 
they describe both the XMMl and XMM2 spectra much 
better than the above common single-component mod- 
els, with the x 2 values decreased by more than 140 



for the total degrees of freedom of about 400 of both 
the XMMl and XMM2 spectra. As a way to model 
the hard component self-consistently, we also fitted the 
spectra wit h SIMPL(MCD) an d SIMPL(BB). SIMPL (in 
XSPEC12; ISteiner et all [20091 is an empirical convolu- 
tion model of Comptonization in which a fraction (/sc) 
of the input seed photons are converted into a power law 
parametrized by an index (Psimpl)- We assume that all 
the scattered photons are up-scattered in energy in this 
model. 

The best-fitting values of the column density are con- 
sistent between XMMl and XMM2, with Npi^.e^'j) 
and (8.6±0.6) x 10 20 cm -2 , respectively, using the model 
MCD+PL. Thus, we chose to fit both spectra with a 
common value of Nr. The final results are given in 
Table [3J The best-fitting values of TVjj arc slightly 
higher than the Galactic value of 6.1 xlO 20 cm~ 2 
from the Leiden/Arg entine/Bonn Survey of Galactic HI 
(|Kalberla et al.ll2005f h probably indicating a small intrin- 
sic absorption. 
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Fig. 4. — The MCD flux versus £;Tmcd using the models 
MCD+PL and SIMPL(MCD). The solid, dotted, and dashed lines 
plot the L oc T 4 , L oc T 3 , and L oc T 5 relations, as a guidance. 



For the model MCD+PL, we show the unfolded spec- 
tra and residuals in Figure [3J In this model, the spectra 
are dominated by the MCD component at energies be- 
low 1 keV. The fraction of the MCD component is about 
84.5% and 96.4% for XMM1 and XMM2, respectively 
(the 0.2-10 keV unabsorbed flux; Table [3]). The 0.2-10 
keV flux increases from XMM1 to XMM2 by a factor of 
8.8 (absorbed) or 5.8 (unabsorbed). We also estimate 
the luminosity, using the bolometric flux of each spec- 
tral component. The disk inclination is uncertain, and 
we assume it to be 60°. The PL component diverges at 
low energies, and we integrate its flux above 0.2 keV. 
We obtain luminosities of 1.70 and 6.38xl0 43 erg s _1 for 
XMM1 and XMM2, respectively (Table©. For compar- 
ison, the corresponding 0.2-10 keV luminosities are 0.47 
and 2.70xl0 43 erg s _1 , respectively. 

We plot the MCD bolometric flux versus its temper- 
ature at the inner disk radius ^Tmcd in Figure [4] (the 
upper panel) . We can see that the evolution of the MCD 
luminosity is consistent with the L oc T 4 track (the solid 
line), which implies a constant inner disk radius with the 
change in luminosity. We note that this is based on the 



only two observations available. The disk temperature 
is relatively low, only /cTmcd = 65.8 and 93.1 eV for 
XMM1 and XMM2, respectively. The PL component is 
weak, and its parameter values have relatively large un- 
certainties. Its index is consistent between XMM1 and 
XMM2 and is relatively steep, with Tpl about 3.5. Forc- 
ing XMM1 and XMM2 to have the same value of r PL in 
the fit, we see a change of the PL normalization iVpL by 
a factor of 2.6 (4.5 a). 

The model SIMPL(MCD) gives results very similar to 
the model MCD+PL (Table |1 and Figure S]), in terms of 
the MCD temperature, the thermal fraction, etc. It in- 
fers that only about /sc=3% and 1% of the thermal disk 
emission is Comptonized to the hard emission in XMM1 
and XMM2, respectively. This model, with a natural cut- 
off at low energies for the hard component, infers lumi- 
nosities similar to those of the model MCD+PL obtained 
by integrating the PL flux down to 0.2 keV (Tabled). 

The spectra can also be fitted almost equally well us- 
ing the models BB+PL and SIMPL(BB) (Table©. The 
BB component dominates in both XMM1 and XMM2, 
contributing >80% of the 0.2-10 keV flux, similar to 
the MCD component in the models MCD+PL and 
SIMPL(MCD). Its effective temperature fcT BB is also 
low, about 58 and 78 eV for XMM1 and XMM2, re- 
spectively. The MCD and BB model s have very sim- 
ilar s pectral shapes at high energies (Mak ishima et all 
1986), but their differences become large at low en- 
ergies. We estimate their differences in the UV. We 
have measurements from two UV filters, i.e., UVW1 
and UVM2, from XMM1. The flux densities of the 
MCD component in the model MCD+PL from XMM1 
are (2.15+0.69) and (3.68+1.18) x 10" 17 erg s" 1 cm' 2 
A' 1 at the effective wavelengths of UVW1 (2910 A) and 
UVM2 (2310 A), respectively. The corresponding val- 
ues for the BB component in the model BB+PL from 
XMM1 are (0.40+0.16) and (1.00+0.40) xlO" 19 erg s" 1 
cm' 2 A' 1 , respectively. The corresponding flux densities 
measured with UVW1 and UVW2 are (3.00+0.47) and 
(5.08+1.45) xl0~ 16 erg s _1 cm~ 2 A -1 , respectively, after 
the Galactic dust extinction c orrection using a red den- 
ing value of £ (B -v) = 0.098 (jSchlegel et al.lll998[) and 
assuming a spectral shape of a MCD model at low fre- 
quencies (i.e., a power law with a photon index of 2/3). 
We see that the UV flux from the OM detection is much 
higher than the BB flux in the UV. It is closer to the 
MCD flux in the UV, but still there is about an order 
of magnitude difference, which will be discussed in Sec- 
tionSJ The flux of the PL component in the UV is hard to 
assess as this model is too steep and diverges at low ener- 
gies (more than two orders of magnitude higher than that 
measured by the OM) and must be cut off below some en- 
ergy. The models SIMPL(MCD) and SIMPL(BB) show 
no such problem, and their fluxes in the UV, from the 
whole model or only from the thermal components, are 
very close to the MCD and BB fluxes in the UV obtained 
above. 

3.3. Fast X-ray Variability 

The left panels of Figure [5] show the pn 250 s 
background-subtracted light curves. The background is 
at a level of about 3% and 2% for XMM1 and XMM2, 
respectively. The variations of the source count rate can 
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5. — Left panels: The pn 250 s 0.2-2.0 keV background-subtracted light curves. Right panels: The PDS of the pn 0.2-2.0 keV light 
curves binned at 199.1 ms for XMM1 and 73.4 ms for XMM2. The black solid constant line is the average PDS above 0.1 Hz, representing 
the Poisson level, and the dotted line is the best-fitting model of a PL plus a constant. 
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TABLE 4 

The fit results of the PDS using a PL plus a constant. 



Obs 


Tpl 


«pl (10- 5 ) 


C P 




rms(%) 


XMM1 
XMM2 


1 - '^-0.57 

1.65±0.19 


3 40+ 27 ™ 

„ ,-+20.42 
S ' 17 -6.75 


25.59±0.17 
2.77±0.01 


1.07(75) 
1.06(82) 


21.2±5.3 
21.1±1.9 



Note. — The column rms refers to 0.0001-0.01 Hz fractional 
rms after subtracting the Poisson level. Wpl is the PL normal- 
ization at 1 Hz. All errors are at a 90%-confidencc level, except 
for the rms, whose 1-cr errors are given. 



be clearly seen for both observations. For XMM2, which 
has higher count rates, we see that the source varies by 
a factor of ^4 within 5 ks. 

The right panels of Figure [5] show the pn PDS of 
XMM1 and XMM2. Both PDS are flat at frequencies 
above 0.1 Hz, representing the Poisson level. Their aver- 
ages above 0.1 Hz weighted by the errors differ from the 
expected Poisson noise values by <0.3%. Below 0.01 Hz, 
both PDS show a clear deviation from the Poisson level 
(the black solid line). We successfully fitted both PDS 
with a power-law (PL) plus a constant Cp, accounting 
for the Poisson level, and the results are given in Table QJ 
The PL index Tpl is about 1.7 for both observations. We 
evaluate the fractional rms within 0.0001 to 0.01 Hz af- 
ter subtracting the Poisson level and obtain a value of 
(21.2±5.3)% and (21.1±1.9)% for XMM1 and XMM2, 
respectively. The above results indicate that the source 
shows similar variability in both XMM1 and XMM2. 

The large fast variability seen above should be due to 
the thermal component from the models in Table as 
the hard component is very weak (e.g., <4% in XMM2, 
in terms of the 0.2-2.0 keV pn count rates). To see the 



cause of the variability, we extracted the high- and low- 
state spectra from the bright observation XMM2, corre- 
sponding to intervals with the pn count rate higher or 
lower than 0.8 counts s _1 , respectively (Figure [5j). We 
fitted these two spectra simultaneously using the model 
SIMPL(MCD) with a common value of N K . We find 
that the MCD temperature in the low state is smaller 
than in the high state by 10 eV at a 4.4-er confidence 
level. Their MCD normalizations are consistent with be- 
ing the same within the error bars. The parameters of 
the hard component have large uncertainties, making it 
hard to constrain any trend. Based on the MCD model, 
the above results provide one explanation of our large 
fast variability: it is caused by the fast variations in the 
mass accretion rate, with the disk truncated at a con- 
stant radius. 

3.4. Comparison with ROSAT and Swift Observations 

We obtain count rates of <1.9xl0 -3 and 1.3 +24 xl0~ 3 
counts s" 1 for the ROSATPSPC observation in 1992 Oc- 
tober and the Swift XRT observation in 2011 February, 
respectively. The 3-cr upper bounds are given above, and 
the lower bounds are zero. Using the response matrices 
corresponding to the source extraction regions used and 
the fit results of the model SIMPL(MCD) for XMM1, 
these count rates correspond to the 0.2-10 keV absorbed 
fluxes of <0.3 and 0.6 +10 xl0 -13 erg s _1 cm~ 2 , respec- 
tively, and the 0.2-10 keV unabsorbed fluxes of <2.4 
and 4.5 +83 xl0~ 13 erg s _1 cm~ 2 , respectively. Thus, 
the source varied by a factor of >64 and >43 between 
the R OS A T pointed observation in 1992 and XMM2, us- 
ing the 0.2-10 keV absorbed and unabsorbed fluxes, re- 
spectively. The corresponding variation factors between 
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XMM2 and the Swift observation in 2011 are >12 and 
>8, respectively. 

4. DISCUSSION 
4.1. The Nature and Implication of the Soft Component 

The remarkable features of 2XMMi J184725. 1-631724 
are the extreme softness of its X-ray spectra and the large 
variability. Understanding the nature and implication of 
the soft component which dominates the X-ray spectra 
will help to pin down the nature of the source. As it is in 
the direction coincident with the center of the galaxy IC 
4765-f01-1504, we first assume that its X-ray emission is 
associated with the SMBH in this galaxy. This can be 
due to either a tidal disruption event or an AGN. 

The fits with the models MCD+PL and SIMPL(MCD) 
above assume that there was a thermal disk emission, 
whose luminosity fraction was inferred to be very high, 
^90%. The evolution of the MCD luminosity is con- 
sistent with L cx T 4 , though only two observations are 
available. These properties are very similar to the ther- 
mal s tate of BH X-ray binaries (|Remillard fc M cClintock 
2006), in which the accretion disk is believed to be trun- 
cated at the innermost stable circular orbit (ISCO). Thus 
we assume that the disk is also truncated at the ISCO 
during these two observations and roughly estimate the 
BH mass from the MCD normalization Amcd- Using a 
distance of 143.9 Mpc and assuming a disk inclination of 
60°, we infer the BH mass to be ^3xl0 5 M©, neglecting 
factors such as the spin and the hardening effect. If we 
replace the MCD model with the more realistic accretion 
disk model ar ound a Kerr black hole kerrbb (in XSPEC; 
iLi et al.l 120051 ) and explore the parameter spaces of the 
disk inclination 0-75°, the spin parameter a* 0-1, and 
the hardening factor 1-1.7, we obtain a range of the BH 
mass of (0.06-3.81) xlO 6 M Q . Assuming the BH mass 
to be 5xl0 5 M , the source would be at about 0.3 and 
1.0 Eddington luminosity in XMM1 and XMM2, respec- 
tively, commo n values seen in t he thermal state of BH 
X-ray binaries ([Done et al.l l2007). We note that the light 
crossing time of the inner accretion disk around a BH 
with this mass is ~50 s, about the timescale on which 
the source begins to show strong variability (Figure [5]). 
The above mass estimate is consistent with that using 
the BH mass vs. bulge luminosity relations (Section l3.ip . 
considering the large uncertainties of both methods. 

Som e AGN can be very soft, showing strong soft ex- 
cesses (jPuchnarewicz et al.l ll992). The soft excess refers 
to the excess of emission below ^2 keV with respect to 
the extrapolation of the power-law fit of the continuum 
above 2 keV and is com monly observed in type-I AGN 
([Turner fe Pounds! 1 1989f ). The strongest soft X-ray ex- 
cesses and variability are found in Narrow Line Seyfert 1 
galaxies (NLSls; e .R.,EQ0m_et^\mB,\LdK^\mM^ 
iGrupe et al1l2010l) . The nature of the soft excess is still 
unclear. The above thermal disk mode l is one of the 
sever al competing models invoked (e.g., iWalter fe Fink! 
1993). The problem with this explanation is that the 
soft excesses from a sample of AGN with a large range 
of mass and luminosity have characteristic temperatures 
spannin g a narrow range (M3.1-0.2 keV), which is hard t o 
explain flGierlinski fe Pond [200l iCrummv et all [2006h . 
The narrow range of the characteristic temperatures of 
the soft excesses finds a natural explanation if they are 



due to atomic processes. The re are two main scenarios , 
i.e., absorption and reflection. iGierlinski fe Done! ([20041 ) 
proposed the soft excess as an artifact of strong, rela- 
tivistically smeared, partially ionized absorption. We 
test this model using the swindl model from XSPEC. 
We use the model swindl (PL), with Tpl required to be 
<3.5 to make sure that it is not the steep PL describ- 
ing the soft excess. We obtain the values of xK 11 ) to be 
1.32(114) and 1.31(283) for XMM1 and XMM2, respec- 
tively. Both observations require strong velocities (>0.5 
and 0.29±0.02 speed of light for XMM1 and XMM2, re- 
spectively) to smear the absorption/emission lines, which 
is har d to achieve from a radiatively driven accretio n disk 
wind dSchurch fe Pond [20071 ISchurch et al.l 12009ft . The 
difference between XMM1 and XMM2 is mainly due to 
different absorption column densities and smearing ve- 
locities, requiring dramatic changes in the absorber. 

In the reflection model, a series of soft X-ray emis- 
sion lines below 2 keV, if strongly relativistically blurred, 
can p roduce the sm o oth s oft excess feature. We 
follow ICrummv et al.l ([2006D to use the model kd- 
blur(PL+reflionx), where kdblur is the relativistic con- 
volution and reflionx is a table model of the ionized 
reflection (see their references therein). In the fits we 
force XMMl and XMM2 to have common values of 
Ah, the inclination, and the Fe abundance. We obtain 
X^)=0.95(lll) and 1.34(281) and reflection flux frac- 
tions ~0.94 and 1 (0.2-10 keV, unabsorbed) for XMMl 
and XMM2, respectively. The fits require a low in- 
clination (17.1±7.6°), a high Fe abundance (8.37t°;™ 
solar value), and a steep illuminating power-law spec- 
trum (reaching the upper index limit of 3.3 allowed in 
the model). The disk emissivity index is different be- 
tween these two observations (5. 37^27 and 8.99±0.89 
for XMMl and XMM2 respectively), implying a very dif- 
ferent disk structure. This model infers a highly spinning 
BH, with the inner disk radius at 2.38±0.38 and 1.87t°;ig 
gravitational radii for XMMl a nd XMM2 respectively . 
Compared with the results from ICrummv et al.l ((2006) , 
the above values are extreme but not unique. 

We see that both the absorption and reflection mod- 
els invoke extreme environments to explain our source. 
The former requires absorbers at very high velocities and 
varying dramatically. The latter requires very steep illu- 
minating spectra and very different disk emissivities be- 
tween XMMl and XMM2. In comparison, the thermal 
disk emission explanation for the soft component in our 
source is more reasonable. In this model, the difference 
between XMMl and XMM2 is simply due to the change 
in the accretion rate. We note that the inferred inner 
disk temperatures are lower than the typical character- 
istic temperatures of the soft excesses in AGN, which 
might indicate that the soft component of our source has 
a different nature from soft excesses in AGN. 

Our spectral fits did not combine both the UV and 
X-ray spectra. There is about one order of magnitude 
difference (about 9a) between the UV fluxes of the ther- 
mal disk inferred from the soft X-ray spectral fits and 
the OM measurements (Section [3]) . This can be due to 
several factors. First, we have only used the simple MCP 
model, while a more realistic disk spectral model is prob- 
ably needed to fit broadband data (about three decades 
in frequency here; lLaor fe N etzcr 198!| IRoss et aLlll992t 




2007 



2009 
Time (yr) 



201 1 



Fig. 6. — The long-term luminosity curve inferred from X-ray 
spectral fits. Arrows represent 3-<r upper bounds. Note that the 
ROSAT observation was made in 1992 October. The solid curve 
is Luminosity = 3.45 X 10 43 [(Time - 2006.60yr)/(lyr)]" 5 '' 3 erg 
s — x , and the dotted curve is Luminosity = 7.57 X 10 43 [(Time — 
2006.18yr)/(lyr)]- 5 / 3 erg s" 1 . 



iVasudevan & Fabian 2009). Secondly, the starlight or 
hot gas emission might be significant in the UV. This is 
supported by the little variability of the UV emission be- 
tween XMM1 and the Swift observation in 2011 Febru- 
ary. Finally, some UV emission might come from the 
reprocessing of the X-ray emission in the outer disk or 
surrounding gas. 

4.2. The Tidal Disruption Event Explanation 

We show above that the X-ray emission in XMM1 and 
XMM2 can be best explained as coming from a ther- 
mal disk around a SMBH with a mass of — 10 5 -10 6 M G . 
The transient nature of the source makes it a great tidal 
disruption event candidate. This is further supported 
by the extreme softness of its X-ray spectra and the in- 
activity of the nucleus of the candidate host galaxy IC 
4765-f01-1504. The inactivity of IC 4765-f01-1504, i.e., 
not an AGN, is indicated by its lack of significant optical 
emission lines (Figure [2]). As our new optical spectrum 
was made four years after the flare, this also implies no 
detection of optical emission lines due to the flare at this 
stage. The 2MASS IR colors (Table [2} also put this 
galaxy in the region occup ied by inactive galaxie s in the 
color-color diagram fe.g- IHvland fc Alien! [1981) . The 
estimated BH mass allows the tidal disruption of a solar- 
type star to be observable (jLidskii fc Ozernolll979l : lReesl 
1988). The luminosity reached 6.4xl0 43 erg s " 1 , which 
is ab out the average seen in other candidates (|Komossal 
12001 lEsquej et all 120071 iGezari et all[2009l ). However, 
no previous candidates had soft X-ray spectra with such 
high quality during the peak of the flare to allow for the 
detailed spectral studies here. 

One main feature of tidal disruption events is the tem- 
poral evolution. Such events are predicted to rise on 
timescales of s everal months and decay on timescales o f 
months/years (jLidskii fe OzernoJll979HReeslll988lll990D . 
The decay approximately follows L oc (t — tn) -5 ^ 3 , where 
£d is the tidal disruption time, for most candidates 



(iMaksvm et al.l 120101: lEscmei et al.l l2008t lHalpern et all 
120041: IK omossa fc Badd Ti999h consistent with the the- 
ory. XMM1 is fainter than XMM2, and they are sepa- 
rated by 211 days. If our source is due to tidal disrup- 
tion, XMM1 should probably be in the rising phase, and 
XMM2 in the decay phase. We plot the luminosity curve 
in Figure [6l Two decay curves following L oc i~ 5 / 3 are 
also plotted, with XMM2 assumed to be in the decay. 
The dotted and solid curves assume £d to be at half a 
year and one month before XMM1, respectively. These 
assumptions are reasonable, as the minimum period for 
the material to return to the SMBH (so as to accrete) 
after disruption is on the order of a month (Reesl fT988c 
lEvans fc Kocha nek 1989]). These curves predict that the 
source luminosity in the Swift observation in 2011 Febru- 
ary should be a factor of about 20 less than XMM2. Our 
detection limit is consistent with this. In the above, we 
concentrate only on the decay, as the connection from 
the rise to the decay is uncertain. The rise could prob- 
ably take months, and there could be a short period of 
several months when the luminosity i s maintained near 
the Eddington limit (|Reesi ri988. 1993). Future sensitive 
X-ray observations in the decay are needed to constrain 
the long-term evolution. 

We see that our source showed large fast X-ray variabil- 
ity, with an rms of -21% (Section [331) in botn XMM- 
Newton observations. If they are really dominated by 
the thermal disk emission, the above rms value is high, 
compared with the values of <5% typically seen in the 
thermal state of neutron- star (NS) or BH low-mass X-ray 
binaries (iLin et al.l 120071: iRemillard fc McClintockl 120061 : 
iMcClintock fc Remillardl 120061 ) . We note that our rms 
was integrated over two decades in frequency, the same 
as the above studies. Besides, in low-mass X-ray binaries, 
the power of the thermal disk typically scales w ith the 
frequency as v~ - 1 (|McClintock fc Remillardl [2006h . while 
the power for our source is steeper, scaling approximately 
as v~ . Large fast variability and steep power can be 
seen sometimes in the thermal state of X -ray binaries, 
such as the flaring branch of bright NSs (jHoman et al.l 
2002), which was ascribed t o som e local instability in 
the inner disk bv ILin et al.l (|2009l ). For the case of our 
source, we have shown in Section 13.31 that its fast vari- 
ability can be explained as due to fast variations in the 
mass accretion rate, which, in the context of a tidal dis- 
ruption event, could reasonably be ascribed to shocks 
during dr astic comp ression and distortion of the stellar 
material (|Reeslll988l ). This does not occur in low-mass 
X-ray binaries, in which the mass is transferred through 
the Roche lobe. Thus we speculate that the large fast 
variability and the steep power in the thermal state is 
intrinsic to tidal disruption events. Among the previ- 
ous tidal disruption event candidates with thermal X- 
ray peaks, the peak of NGC 5905 was the best observed. 
It showed an increase of a factor of —3 du ring the peak, 
but it was over four days (|Bade et al.lll996l ). and the data 
quality was not high enough to investigate the variability 
on much shorter timescales. 

4.3. Comparison with AGN 

2XMMi J184725. 1-631724 is unlikely to be an AGN 
from the lack of bright optical emission lines. In the fol- 
lowing we briefly compare the X-ray properties of our 
source with ultrasoft AGN to show their similarities and 
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differences. NLS1 galaxies typically show the steepest 
soft-X-ray spectra amon g AGN and have a typical pho- 
ton i ndex of around 3 ([Boiler et al.l 119961 : iGrupe et al.l 
120101 ). The simple fits of the XMM1 and XMM2 spectra 
using the PL model indicate an extreme softness of our 
source that is not seen in NLSls. 

We estimate the optical-to- X-ray spectral slope a ox 
using the flux densities in the rest frame of 2500 
A and 2 keV (jTananbaum et al.l 119791 ) . We obtain 
Qiox— 1-76±0.03. The luminosity density in the rest frame 
of 2500 Ais/ 25 oo=(l-48±0.24)xl0 28 ergs- 1 Hz" 1 . With 
this value of /2500 , the value of a x should b e 5,1-3, based 
on th e sample of 92 Seyfert 1 galaxies in IGrupe et al.l 
([20101 their Figure 15). In fact 2XMMi J184725. 1-631724 
has about the highest value of a ox and the lowest value 
of ^2500 and is an outlier, compared with their sample. 
We note, however, the possible large uncertainty of a ox 
for our source, whose star light contamination might be 
large, as discussed above. 

There have been several ultrasoft AGN claimed in the 
literature. The NLS1 galaxy WPVS 007 has the softest 
X-ray spectrum among AGN detected during the ROSAT 
All-Sky Surve y, with r PT ~8 if fitted with a PL or fcT~20 
eV with a BB (|Grupe et al.l[l995l). It can be exp lained as 
emission from the inner disk (jGrupe et ailll995l ) , but the 
quality of the data and the lack of simultaneous observa- 
tions above 2.4 keV could mak e its soft excess due to the 
presence of a warm absorber ([Grupe et al.l 1200 8b). The 
narrow-line quasar PHL 1092 has XMM-Newton obser- 
vations, and the X-ray spectra are stee p (rpT,~4-5), but 
the PL fit is n ot goo d when it is bright (jGailo et alll2004t 
iMiniutti et~ail2009D . Its bright spectrum in 2003 was fit- 
ted with a model of MCD+ PL plus an absorp tion line 
and a reflection component bv lGallo et al.l (|2004f ) : its soft 
excess was ascribed to the MCD component (fcT=114±4 
eV), as in our study. The NLS1 galaxy 1H 0707-495 
also shows an intense soft exce ss. Its X-ray spect ra from 
XMM-Newton show r PL ~3.8 ([Boiler et all 120021 ) . They 
were fitted with a model of BB+ PL plus a reflection com- 
ponent by iFabian et al.l (|2009l) . From this model, both 
the thermal disk emission and the reflection contribute to 
the soft excess, but the latter dominates. We note that 
the above ultrasoft AGN show large short-term variabil- 
ity factors of a few and/or long-term variability factors of 
a few hundred, similar to our source. Thus the variabil- 
ity of our source is not extreme compared with NLSls, 
but the softness of its X-ray spectra is hardly challenged. 
We note that the soft component in our source cools in 
low states, which is not generally seen in NLSls (e.g., 
IMiniutti et"aTll2009D . 

4.4. Alternative Explanations 

We explore the possibility that 2XMMi J184725.1- 
631724 is not associated with the SMBH in IC 4765-f01- 
1504. If it is an ultraluminous X-ray source (ULX) in 
this galaxy, it would have a luminosity about one order 
of magnitude brighter than the brightest ULX reported 
thus far, i.e., HLX-1 (jFarrell et al.l 12009) . However, as 
derived above, both the spectral fits and the variability 
argument imply that this source most probably has a 
mass of ~10 5 -10 6 M , and it should be within 0.5 kpc 
(3-cr error) of the galaxy center. Being so massive and 
so close to the galaxy center, it seems unlikely to be a 



source other than the central SMBH. Thus we deem that 
our source is not a ULX in IC 4765-f01-1504. 

The other possibility is that 2XMMi J184725.1- 
631724 is a foreground Ga lactic source. Following 
lHaakonsen fc Rutle dgc (|2009fl and using the 0.2-10 keV 
absorbed flux of 1.92X10" 12 erg s" 1 cm" 2 (Table©, 
we have the X-ray to IR flux ratio fx/ fi ~ 5 for 
our source, making it unlikely to be a coronally ac- 
tive star. These stars gene rally have fx/ fi < 0.03 
(jHaakonsen fc Rutledge1l2009D . 

The softness of the source with characteristic tempera- 
tures of a few tens of e V makes it similar to the super-soft 
X-ray source s (SSS; IKahabka fc van den Heuvell [2006; 
IGreinerl [20 00) . This class of objects have BB temper- 
atures in the range 20-100 eV, which are about two or- 
ders of magnitude lower than X-ray binaries contain- 
ing an accreting NS or BH. There are various types 
of SSS. One main class is coolin g white dwarfs (WDs; 
IKahabka fc van den Heuvell [2001 . 2XMMi J184725.1- 
631724 brightened by a factor of >64, thus ruling out 
this hypothesis. A large fraction of SSS can be in- 
terpreted as nuclear burning of the hydrogen-rich mat- 
ter on the surface of a white dwarf that accretes 
matter from the companion, such as in the so-called 
close binary super-soft sources and super-soft novae 
([Kahabka fc van den Heuvel [2001 IGreinerl [2000). It 
seems unlikely that our source is such based on the fol- 
lowing considerations. In our Galaxy, there are only 
about a dozen such source s detected sinc e the Einstein 
Observatory observations ([Greinerl H000) , indicating a 
very low density of such objects in the sky or a very 
low life duty cycle. The probability for any of them ly- 
ing in the direction of the center of a galaxy is simply 
negligible. Besides, these objects are typically observed 
at luminosities of ~10 36 -10 38 ergs s _1 , while our source 
has much lower luminosities, <10 34 ergs s _1 , if it is 5 kpc 
awa y. These objec ts are mostly found at distances of <5 
kpc (jGreinerll2000T) . 



5. CONCLUSION 

2XMMi J184725. 1-631724 is an ultrasoft X-ray tran- 
sient source with characteristic temperatures of a few 
tens of eV. It was bright in two XMM-Newton observa- 
tions in 2006-2007, but was not detected in a ROSAT 
pointed observation in 1992, implying a variation factor 
of >64 in the 0.2-10 keV absorbed flux. It was unde- 
tected again in a Swift observation in 2011 February, im- 
plying a flux decrease by a factor of >12. It lies toward 
the center of the galaxy IC 4765-f01-1504 at a redshift 
of 0.0353. No bright optical emission lines were detected 
from this galaxy, making this source a good tidal disrup- 
tion event candidate. The fits to the two XMM-Newton 
spectra using a thermal disk plus a weak hard compo- 
nent indicate that the accretion disk luminosity appears 
to follow the L cx T 4 relation and that the BH mass is 
around 10 5 -10 6 M©. The source showed large fast vari- 
ability in both XMM-Newton observations, which can be 
explained as due to fast variations in the mass accretion 
rate. To further check whether this is a tidal disruption 
event, future long-term X-ray monitoring is necessary to 
see whether it follows the decay expected for a tidal dis- 
ruption event. 
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